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Abstract 

A variational analysis is presented for the generalized spiked harmonic oscillator 
Hamiltonian operator — ^ + Bx 2 + £ + ^ , where a is a real positive parameter. The 
formalism makes use of a basis provided by exact solutions of Schrodinger's equation 
for the Gol'dman and Krivchenkov Hamiltonian, and the corresponding matrix ele- 
ments that were previously found. For all the discrete eigenvalues the method provides 
bounds which improve as the dimension D of the basis set is increased. Extension to 
the TV -dimensional case in arbitrary angular-momentum subspaces is also presented. 
By minimizing over the free parameter A, we are able to reduce substantially the 
number of basis functions needed for a given accuracy. 
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I. Introduction 

Since the fascinating work of Harrell [1] on the ground-state energy of the singular 
Hamiltonian H = H Q + XV = -d 2 /dx 2 + x 2 + X/x a , x e [0, oo], a > 0, known as the spiked 
harmonic oscillator Hamiltonian, the volume of research in this field has grown rapidly. 
This is not only because of the important applications of singular Hamiltonians to a 
wide variety of problems in chemical, nuclear and particle physics, but also because of 
its intrinsically interesting properties from the point view of mathematical physics [2- 
5]. Most of these work [6-15], however, has focused on studying the spiked harmonic 
oscillator Hamiltonian in one spatial dimension since the interesting Klauder phe- 
nomenon [2-5] associated with H does not occur in higher dimensions. Klauder [2-5] 
has shown that, for sufficiently singular potentials, V cannot be turned off smoothly 
in the Hamiltonian H to restore the free Hamiltonian H . Aguilera-Navarro et al [6] 
employed variational and perturbative schemes to solve the spiked harmonic oscillator 
problem for the ground state energy. In their variational analysis of the Hamiltonian 
H , Aguilera-Navarro et al employed the basis set of harmonic oscillator eigenfunc- 
tions normalized in the interval [0,oo] , i.e. the set of Hermite functions generated by 
the non-singular harmonic-oscillator potential x 2 . 

Recently, we have obtained closed-form expressions [16] for the singular-potential 
integrals < m\x~ a \n > using the Gol'dman and Krivchenkov eigenfunctions [17] for 
the singular Hamiltonian 

H Q = --^ + Bx 2 + 4, B>0,A>0. (1.1) 

ax x 

We present a variational analysis of the generalized spiked harmonic oscillator Hamil- 
tonian 

H = H + \V = -^L + Bx 2 + ^ + ±, (1.2) 

where A and a are positive parameters. To evaluate the matrix elements of x~ a , Hall 
et al [16] used the basis set constructed with the normalized solutions of Schrodinger's 
equation H a tp = Etp , i.e. 

il> n (x) = \n> = C n xM 1+ ' /r & X )e-t VB * 2 iF 1 (-n, 1 + ^VT+4A; VBx 2 ); 

2BH^™r(n + l + lyT+42) ( L3 ) 
n![r(l + ± v / T+43)P ' n - U '^"-" 

where ii<i is the confluent hypergeometric function [18] 

ifi(a, b;z) = }2 ^g^p (o) fc =a(a + l)...(a + k-l)= ^ 
Hall et al found, for a < 27, that the matrix elements <m\x~ a \n> are given by 



<m\x \n> 




T(7 + m) 



n\m\ T(7 + n) 

(1.4) 

r(fc + 7 -f)r(f-A! + n) 



r(fc + 7 )r(f-fc) ' ' 2 



-V1 + 4A, 
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in which each element has a factor which is a polynomial of degree m + n in a . Of 
particular interest are the matrix elements < 0\x~ a \n > , which can be obtained from 
Eq.(1.4), and read as follows: 



<0|z \n>-{l)B y n , r(7 + n) r( 7 )r(f) >»-°>M,--.- 

The special case where A = and B = 1 allows us to recover the matrix elements 
provided by Aguilera-Navarro et al [6] for the operator x~ a in the harmonic oscillator 
representation, supplemented by the Dirichlet boundary condition ^(0) = . Indeed, 
substituting A = and B = 1 , and using the identities 

r(n+|) (2n+l)! r(z + n) , 1W „, , 

after some algebraic simplification, we can easily show that Eq.(1.4) reduces to Eq.(12) 
and Eq.(13), for m = , of Ref. [6]. All these expressions are valid for a < 3 . 

The purpose of this paper is to employ the variational method, with the matrix 
elements (1.4), to solve the generalized singular Schrodinger equation 



d 2 2 A A 
dx^ + Bx+ V 2+ ^ 



ip = Etp, < x < oo. (1.6) 



The present work is a generalization of the variational approach of Aguilera-Navarro 
et al [6] that used harmonic-oscillator functions and was restricted to ground state 
level (B = 1,A = 0). 

The paper is organized as follows. In Sec. II, we outline the variational method 
used to study (1.6), and we also extend the scope to cover the N~ dimension case. 
Some numerical results, and comparisons with the results of Aguilera-Navarro et al, 
are presented in Sec. III. In Sec. IV it is shown that a further optimization over 
the free parameter A reduces substantially the number of basis functions needed to 
compute the eigenvalues to a given accuracy. 

II. The variational method 

The first step in the variational method is to select a suitable complete set of 
basis functions that is adapted to the problem at hand. In the variational analysis 
of the ground state energy of the singular potential V(x) = x 2 + \x~ a , known in the 
literature as the spiked harmonic oscillator potential, Aguilera-Navarro et al employed 
a basis set of harmonic oscillator eigenfunctions normalized on the interval [0, oo] and 
vanishing at x = , i.e. the set of odd Hermite functions generated by the non-singular 
harmonic-oscillator potential x 2 . A more effective basis set for the variational analysis 
of such singular problems is the set of normalized wavefunctions (1.3) because the 
singular characteristics of the potential are naturally built into the wavefunctions. 
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Let ip(x) be a 'trial function' for Hamiltonian H given by (1.2), and let us suppose 
that ip(x) can be expanded in terms of the basis set ip n (x) defined by (1.3). Thus we 
have 

D-l 

^{x) = ^2 a ni } n{x). (2.1) 
n=0 

The problem now is to minimize the eigenenergies of (1.2), with respect to the vari- 
ational parameters a n , n = 0, 1, ...,D — 1 , in the finite dimensional subspace H D 
spanned by the D functions ipo,4>i, ■ ■ ■ ,4>d-i ■ This variational problem is equivalent 
to diagonalizing the Hamailtonian (1.2) in the subspace H D . By increasing the di- 
mension D, we can always improve the results. Thus, we have to evaluate the matrix 
elements of the Hamiltonian (1.2) in the basis (1.3). They can be separated into two 
contributions 

H mn = <m\H\n> = <m\H \n> + \<m\x~ a \n>, m, n = 0, 1, 2, ...,D- 1. (2.2) 

Since H is diagonal in the chosen basis, the first term on the right hand side of (2.2) 
is the exact solutions of the Gol'dman and Krivchenkov Hamiltonian, that is 

<m\H a \n> = V~B(4n + 2 + y/l + 4A)S mn , m, n = 0, 1, 2, ...,£>- 1, (2.3) 

where S mn is the kronecker delta that equal to 1 if m = n and if m ^ m . The 
second term is given by the matrix elements (1.4). Explicit expressions for the first 
fifteen matrix elements of x~ a are given here in the Appendix. 

In order to extended the scope of this analysis to the N -dimensional Schrodinger 
Equation (1.8), we observe first that the A term has the dimensions of kinetic en- 
ergy, such as the term that appears in higher-dimensional systems. We may therefore 
replace A in Eq.(2.2) with 

A^A+{l + ±(N-l))(l+±(N-3)), N>2, (2.4) 

where the unperturbed energy levels (2.3) becomes in this case 



<TO|#ol«> = 2v^2n+l + yyL+(/ + y-l) 2 )<W, N>2. 



Thus, Eq.(2.2) becomes 



H mn = 2sfB\ 2n+l + \A+{l + ^ - l) 2 )<5 mn + A<m|x- Q |n>, m, n = 0, 1, 2, ...,£>- 1, (2.6) 



V V 2 

where the matrix elements <m\x~ a \n> becomes 



<m\x~ a \n> = {-l) n+m B- 



T{m+l + ^A+(l + N/2-lf) 
n!m! I> + 1 + y/A+(l + N/2-1) 2 ) 

A+(l + N/2-iy)T{'±-k + n) 



y f u k(m y(k-§ + l + VA 
to W T(k + l + VA 



+ (i+N/2-mm-k) 
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We have omitted the case N — 1 because this curious singular problem in one dimen- 
sion has features [19-21] that are not in harmony with our main purpose. 

To recover the results of Aguilera-Navarro et al, we substitute A = , B = 1 , 
N = 3 ,and I = in Eq.(2.6). In other words, we obtain a general variational expression 
(2.2) that treats the solution of Schrodinger's equations (1.8), and also the work of 
Aguilera-Navarro et al in a single formulation. 

III. Some numerical results 

From Eqs. (2.2), (2.3), and (1.4) and the results given in the Appendix, it can be 
readily seen that the first variational approximation (subspace of dimension 1) to the 
ground state eigenvalues of the Hamiltonian (1.2) is 



E a = #oo = <0\H |0> = + 



r(7-f) 
r(7) : 



7 = 1 + -V1 + 4A 



For the case A = and B = 1 , E reduces to 



^(o) — 3 + A 



r(f) 



(3.1) 



(3-2) 



which coincides with the ground-state energy expression for the spiked harmonic 
oscillator obtained by Aguilera-Navarro et al, i.e. Eq.(5.1) in Ref.[6]. 

When D = 2 the diagonalization can also be performed analytically, by means of 
the secular equation, and we obtain 



E± = - 
* 2 



4V^(1 + 7) + \b? (a 2 -2a + 87) ^^ 



A6B + 2XB^a{a-2) 



Th - 2l) A 2 

:> +— Bfa 2 (( a -2) 2 + 16 7 ) 



r(7-f) 



r(7 + i) 



(3.3) 



r(7 + i) 

where E = E_ and E 1 = E + . Again, if we put A = 0, B = 1, and a = 5/2 we obtain 

< 3-a \ 



E ± ='- 
* 2 



10 



53 r( 2 
24 A r(f) 



± 4/I6+ -A 



5 x r(^s 



r(|) 



£(fg) 
r(l) 



(3.4) 



which coincides with Eq.(5.2) in Ref.[6]. 

For higher values of the variational space dimension D we have to use a numerical 
diagonalization procedure such as that of Jacobi [22] to find the eigenvalues of the 
matrix H given by 









Hqd-1 


\ 






Hn . 


HiD-l 




\H D-10 


Hd-ii 


Hd-id- 


1/ 



(3.5) 



Table I shows the convergence of these results, when a = 2.5 , for the ground state 
energy (i.e. A = 0, B = 1 ) of the spiked harmonic oscillator, and selected values of A . 
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For this special case with A = and a = 5/2 , our results confirm those of Ref. [6] (their 
Table I), except for minor rounding errors: it appears that the results of Aguilera- 
Navarro et al have been truncated rather than rounded. It is important to mention 
here that, although the matrix element Xg 3 Q , reported in Ref. [6], contains some errors 
in the coefficients for the exponent a, the results in Table I of Ref. [6] are correct. 
Indeed, x^ should read: 



x- Q = -j-^-i( a 6 - 6a 5 + 106a 4 - 384a 3 + 2080a 2 - 3408a + 5040) (3.6) 
7.1 (j) 



instead of 



, (a 6 - 6a 5 + 106a 4 - 454a 3 + 1660a 2 - 3968a + 5040) (3.7) 

7.1 (2) 

as quoted by Aguilera-Navarro et al [6]. The results quoted in their Table I must 
therefore have been calculated with the correct formula. We have also corrected in 
the Appendix here some errors in the coefficients of the general matrix element x^ 
presented in Ref. [16]. 

In Table (II) here we report the upper bounds E nt = Eq obtained for the Hamil- 
tonian H = -■£? + x 2 + ^ in spatial dimensions N = 2 to N = 10. Now the value 
of A depends on the angular momentum I and the number of spatial dimension N 
as determined by Eq.(2.4). As we mentioned above, the method we have discussed in 
Sec. II provides bounds on all the eigenvalues in a given angular momentum subspace. 
We present in Table (III) our results for the eigenvalues E nt = E^ for spatial dimen- 
sions N = 2 to ./V = 10; the potential and the corresponding wave functions are shown 
in Fig.(l). For comparison, we have integrated Schrodinger's equation numerically, 
with the Hamiltonian given by (1.2), for different values of exponent a. To avoid 
difficulties caused by the singular character of the potential near the origin, we begin 
the integrations near the potential minimum and integrate in both directions, away 
from the starting point. A similar approach has been described by Diaz et al [23]. 

IV. A further variational refinement 

We now introduce another variational adjustment that will substantially reduce 
the number of the basis elements required to compute the eigenvalues of the Hamil- 
tonian, i.e. 

d 2 (1 + ±(JV-l))(Z + ±(i\r-3)) d2 A 
H = 'd^ + x^ + BX + ^ (41) 

We notice first that the Hamiltonian (4.1) can be written as 

H- d " 1 V + h(N-mi + h(N-3)) 2 , A [X A\ 

H - —fa2 + ~2 + BX + ^2 + I — - ^2 ) < 4 ' 2 > 



x* r \s 



where A is additional variational parameter, different from zero, to be determined 
later. In this case Eq.(2.2) becomes 

H mn = 2VB(2n + 1 + y/A + (I + N/2 - l) 2 )<J mn + \<m\x- a \n> - A<m\ X - 2 \n>, (4.3) 
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where <m\x~ a \n> is given by (2.7) and <m\x~ 2 \n> is obtained by setting a = 2 in 
this formula. Upper bounds to the eigenvalues E% are again provided by finding the 
eigenvalues of the matrix H given by Eq.(3.5), but with the entries depending on A 
according to (4.3). 

For example, when the variational space has dimension D = 1, the lowest eigen- 
value of the problem in TV -dimensions labelled by I is determined by 

ffoo = v^(2 7 - -i-) + Aff"/ 4 T{ 1 pi } , (4.4) 
7-1 1 7 



where 7 = l + ^/A + (; + f - l) 2 . For N = 3,1 = 0, the minimum of H m with respect to 
A for a = 2 occurs at A = A , thus yielding the exact solution of the spiked harmonic 
oscillator E = VB(1 + VI + 4A) . In fact, a 'good' general estimate for the value of A 
is A = A. This rough estimate for A reduces substantially the number of the basis 
function needed to compute the eigenvalues; by minimizing over A , we obtain even 
better upper bounds, as Table (IV) clearly indicates. 

V. Conclusion 

We have generalized the work of Aguilera-Navarro et al to treat the more general 
spiked harmonic oscillator problem (1.2). In particular, we have presented a variation 
method to solve the interesting N -dimensional spiked harmonic oscillator problem 
with the x~ a singular term. The present work is limited by the necessary condition 
a < 3. Some interesting results for a > 3 may be found, for example, in Refs.[24- 26]. 

We hope that our work will encourage further research into the spectra generated 
by this interesting class of singular potentials. It is very clear from our results that 
use of a variational basis which is itself derived from a related soluble singular prob- 
lem leads to very effective approximation methods for more general problems of this 
singular type. The presence of the free parameter A in the class of soluble problems 
allows a further refinement in the upper energy estimates. 

Acknowledgment 

Partial financial support of this work under Grant No. GP3438 from the Natural 
Sciences and Engineering Research Council of Canada is gratefully acknowledged. 



Variational analysis for 



page 8 



References 

[1] E. M. Harrell, Ann. Phys. 105, 379 (1977). 

[2] J. R. Klauder, Acta Phys. Austriaca Suppl. 11, 341 (1973). 

[3] L. C. Detwiler and J. R. Klauder, Phys. Rev. D 11, 1436 (1975). 

[4] B. Simon, J. Functional Anal. 14, 295 (1973). 

[5] B. DeFacio and C. L. Hammer, J. Math. Phys. 15, 1071 (1974). 

[6] V. C. Aguilera-Navarro, G.A. Estevez, and R. Guardiola, J. Math. Phys. 31, 99 
(1990). 

[7] V. C. Aguilera-Navarro and R. Guardiola, J. Math. Phys. 32, 2135 (1991). 

[8] Solano- Torres, G. A. Esteves, F. M. Fernandez, and G. C. Groenenboom, J. Phys. 
A: Math. Gen. 25, 3427 (1992). 

[9] M. F. Flynn, R. Guardiola, and M. Znojil, Czech. J. Phys. 41, 1019 (1993). 

[10] M. Znojil, Proc. Int. Conf. on Hadron Structure 91, Proc. CSFR vol.1, 1 (1992). 

[11] M. Znojil, J. Math. Phys. 34, 4914 (1993). 

[12] E. S. Estevez-Breton and G. A. Estevez-Breton, J. Math. Phys. 34, 437 (1993). 

[13] R. Hall and N. Saad, Can. J. Phys. 73, 493 (1995). 

[14] R. Hall and N. Saad, J. Phys. A: Math. Gen. 31, 963 (1998). 

[15] R. Hall and N. Saad, J. Phys. A: Math. Gen. 32, 133 (1999). 

[16] R. Hall, N. Saad and A. von Keviczky, J. Math. Phys. 39, 6345-51 (1998). 

[17] I. I. Gol'dman and D. V. Krivchenkov, Problems in Quantum mechanics (Perga- 
mon, London, 1961). 

[18] L. J. Slater, Confluent Hypergeometric Functions (At the University Press, Cam- 
bridge, 1960). 

[19] K. Bhattacharyya and R. K. Pathak, Int. J. Quantum Chem. 59, 219-226 (1996). 
[20] C. J. Ballhausen, Chem. Phys. Lett. 146, 449 (1988). 
[21] P. Senn, Chem. Phys. letters 154, 172 (1989). 

[22] H. W. Press, P.B. Flannery, A. S. Teukolsky and T. W. Vetterling, Numerical 
recipes in Pascal: The art of scientific computing (Cambridge University Press, 
Cambridge, 1989). 

[23] Carlos G. Diaz, Francisco M. Fernandez and Eduardo A. Castro, J. Phys. A: Math. 
Gen. 21, Lll (1988). 



Variational analysis for 



page 9 



[24] E. Papp, Europhys. Lett. 9, 309 (1989). 
[25] M. Znojil, Phys. Lett. A 255, 1 (1999). 
[26] M. Znojil, Phys. Lett. A 259, 220 (1999) 



Variational analysis for 



page 10 



Appendix: Some explicit forms of the matrix elements <m\x ~ a \n> 

We present the first fifteen matrix elements of x~ a that we used to compute the 
variational eigenvalues in the next sections. Some errors in the coefficients of in 
Ref. [16] have been corrected. In terms of the parameter 7 = 1 + \\J\ + 4A , the explicit 
matrix elements, from Eq.(1.4) and Eq.(1.5), are as follows: 

£«m = <m\x- a \n> 

a _ „ f r(-f+ 7 ) 
^oo - Di r( 7 ) 

x m a = -Sf <* r( " 2+7) 



01 - - 2^ r( 7 ) 

a(a+2) r(-f+- 



'-02 



„-a = «(q+2)(«+4) r(-f+ 7 ) 

' 03 2-V3h(7+l)(7+2) r W 

„-a _ a(a+2)(a+4)(q+6) r(-f+ 7 ) 

J ° A 2V4!7(7+l)(7+2)(7+3) T( -^ 

-a. _ n f (a 2 -2a+4 7 ) T(-f + 7 ) 

-11 — i 1 r( 7 +i) 



X 12 — 



-13 



X 00 — 



p°L a(a 2 -2a+8 7 ) r(-f +7) 
2 3 A /2!( 7 +l) r (7+l) 

a(a+2)(a 2 ~2a + 12 7 ) T(-f + 7 ) 
2V3!( 7 +l)(7+2) T( 7 +l) 

p° «(a+2)(«+4)(a 2 -2a + 167) r(-f +7) 
2S x /4!(7+l)( 7 +2)( 7 +3) r (T+l) 

ggL a 4 -4a 3 + (12+167)a 2 -(16+327)a+327(l+7) r(-f +7) 



22 " "" 2 4 v^!2l r( 7 +2) 

-a _ nf a(a 4 -4a 3 + (20+247)a 2 ~(32+487)a+967(l+7)) T(-f +7) 
^ 23 _ 2V2!3!(7+2) r ^+ 2 ) 



2^94 — 



g2. a(a 5 -2a 4 + (20+327)a 3 +8a 2 + (-96+647+1927 2 )a+3847(7 + l)) T(-f +7) 
2« v '2!4!( 1 +2)( 7 +3) r <7+2) 



3:33° = ii ^L ? ^a 6 - 6a 5 + (52 + 36 7 )a 4 - (168 + 144 7 )a 3 + (352 + 720 7 + 288 7 2 )a 2 - (384+ 1152 7 + 
576 7 2 )a + 384 7 (1 + 7 )(2 + 7)) r fc f + ^ } 



B4 



a 7 - 6a 6 + (76 + 48 7 )a 5 - (264 + 192 7 )a 4 + (832 + 1536 7 + 576 7 2 )a 3 - (1 152 + 



27,/3!4!(7+3) 

v 2\„,2 , 1 rco R „,/„, _i_ 1 V„, _i_ o\™ \ r (~f +7) 



2688 7 + H52 7 2 )a 2 + 1536 7 ( 7 + 1)( 7 + 2)a ) 



" 44 2 s ;S 



a 8 - 8a 7 + (136 + 647)a 6 - (704 + 384 7 )a 5 + (3856 + 4480 7 + 1152 7 2 )a 4 - (10880 + 
15360 7 + 4608 7 2 )a 3 + (19200 + 48640 7 + 32256 7 2 + 6144 7 3 )a 2 - (18432 + 67584 7 + 55296 7 2 + 
12288 7 3 )a + 6144 7 ( 7 + 1)( 7 + 2)( 7 + 3)) \!*$ 
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Table(I) The ground-state eigenvalue of Schrodinger's equation Hip = Etp , where 
A = and H = — ^ + x 2 + , obtained by diagonalization of the D x D matrix ele- 
ments with D = 1, 2, 10, 20, 30 . The "exact" values E were obtained by direct numerical 
integration of Schrodinger's equation. 



A 


£7(1) 


£7(2) 


£7(10) 


£7(20) 


£7(30) 


E 


1000 


4094.062 692 


324.897 482 


44.967 048 


44.955 485 


44.955 485 


44.955 485 


100 


412.106 269 


36.802 319 


17.541 891 


17.541 890 


17.541 890 


17.541 890 


10 


43.910 627 


7.951 034 


7.735 637 


7.735 135 


7.735 114 


7.735 111 


5 


23.455 313 


6.304 224 


6.297 319 


6.296 710 


6.296 566 


6.296 473 


1 


7.091 063 


4.688 098 


4.354 248 


4.329 430 


4.323 263 


4.317 312 


0.5 


5.045 531 


4.216 200 


3.919 692 


3.882 149 


3.869 547 


3.848 553 


0.1 


3.409 106 


3.366 867 


3.316 061 


3.302 484 


3.296 024 


3.266 871 


0.05 


3.204 553 


3.193 800 


3.177 840 


3.172 753 


3.170 127 


3.152 420 


0.01 


3.040 911 


3.040 476 


3.039 702 


3.039 409 


3.039 244 


3.036 665 


0.005 


3.020 455 


3.020 346 


3.020 148 


3.020 071 


3.020 027 


3.019 086 


0.001 


3.004 091 


3.004 087 


3.004 079 


3.004 075 


3.004 074 


3.004 014 



Variational analysis for 



page 12 



Table(II) Upper bounds Eq for A = and H = —-^ + x 2 + for dimension 
N = 2 to 10 , obtained by diagonalization of the D x D matrix H , D = 1,2, 10, 20, 30 . 
The "exact" values E were obtained by direct numerical integration of Schrodinger's 
equation. 



N 


£7(1) 


E (2) 


£7(10) 


£(20) 


£(30) 


E 


2 


196.700 853 


9.092 284 


8.485 580 


8.485 399 


8.485 384 


8.485 378 


3 


21.236 010 


8.698 978 


8.564 442 


8.564 364 


8.564 358 


8.564 356 


4 


13.735 043 


8.813 825 


8.795 449 


8.795 440 


8.795 440 


8.795 440 


5 


11.686 537 


9.163 174 


9.163 095 


9.163 094 


9.163 093 


9.163 093 


6 


11.110 897 


9.650 211 


9.646 713 


9.646 702 


9.646 701 


9.646 701 


7 


11.145 653 


10.233 096 


10.225 061 


10.225 046 


10.225 045 


10.225 045 


8 


11.492 447 


10.889 178 


10.879 092 


10.879 078 


10.879 077 


10.879 077 


9 


12.020 404 


11.603 187 


11.592 993 


11.592 982 


11.592 982 


11.592 982 


10 


12.662 990 


12.363 513 


12.354 191 


12.354 183 


12.354 183 


12.354 183 
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Table(III) Upper bounds E 21 for A = and H = -A + x 2 + for dimension N = 2 
to 10, obtained by diagonalization of the D x D matrix H for D = 30. The "exact" 
values E 2 i were obtained by direct numerical integration of Schrodinger's equation. 



N 




E21 


2 


16.543 648 


16.543 629 


3 


16.904 446 


16.904 445 


4 


17.381 709 


17.381 708 


5 


17.955 446 


17.955 444 


6 


18.607 070 


18.607 067 


7 


19.320 693 


19.320 691 


8 


20.083 407 


20.083 406 


9 


20.885 022 


20.885 021 


10 


21.717 608 


21.717 608 



Variational analysis for 



page 14 



Table(IV): Eigenvalues of Schrodinger's equation Hip = Eip , where H = -j^ + x 2 + ^ 
obtained by diagonalization of the n x n matrix H , n = 1,2, ... ,5 and minimizing over 
the parameter A . 



A 




E (2) 


£7(3) 


£7(4) 


£7(5) 


E 


1000 


44.959 423 


44.955 722 


44.955 562 


44.955 559 


44.955 517 


44.955 485 


100 


17.546 305 


17.542 630 


17.542 082 


17.542 066 


17.542 040 


17.541 890 


10 


7.40 873 


7.736 864 


7.735 869 


7.735 645 


7.735 596 


7.735 111 


5 


6.302 942 


6.298 821 


6.297 638 


6.297 281 


6.297 145 


6.296 473 


1 


3.325 682 


4.321 615 


4.320 076 


4.319 376 


4.318 963 


4.317 312 


0.5 


3.857 330 


3.853 611 


3.852 085 


3.851 313 


3.850 823 


3.848 553 


0.1 


3.273 542 


3.271 566 


3.270 632 


3.270 082 


3.269 700 


3.266 871 


0.05 


3.157 126 


3.155 956 


3.155 378 


3.155 022 


3.154 768 


3.152 420 


0.01 


3.037 845 


3.037 674 


3.037 581 


3.037520 


3.037 474 


3.036 665 


0.005 


3.019 610 


3.019 553 


3.019 522 


3.019 500 


3.019 484 


3.019 086 


0.001 


3.004 053 


3.004 050 


3.004 049 


3.004 048 


3.004 047 


3.004 014 
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Figure (1) The potential y = V(r), and the nine wave functions corresponding to 
the eigenvalues E nl = E 2 i for spatial dimensions N = 2 to 10. 



